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ABSTRACT 

We present numerical hydrodynamical simulations of the formation, evolution and gravitational 
collapse of isothermal molecular cloud cores induced by turbulent compressions in spherical geometry. 
A compressive wave is set up in a constant sub- Jeans density distribution of radius r = 1 pc. As 
the wave travels through the simulation grid, a shock-bounded spherical shell is formed. The inner 
shock of this shell reaches and bounces off the center, leaving behind a central core with an initially 
almost uniform density distribution, surrounded by an envelope consisting of the material in the shock- 
bounded shell, with a power-law density profile that at late times approaches a logarithmic slope of —2 
even in non-collapsing cases. The central core and the envelope are separated by a mild shock. The 
resulting density structure resembles a quiescent core of radius < 0.1 pc, with a Bonnor-Ebert-like 
(BE-likc) profile, although it has significant dynamical differences: it is initially non-self-gravitating 
and confined by the ram pressure of the infalling material, and consequently, growing continuously in 
mass and size. With the appropriate parameters, the core mass eventually reaches an effective Jeans 
mass, at which time the core begins to collapse, until finally a singularity is formed. Thus, there is 
necessarily a time delay between the appearance of the core and the onset of its collapse, but this is not 
due to the dissipation of its internal turbulence as it is often believed, but rather to the time necessary 
for it to reach its Jeans mass. These results suggest that pre-stellar cores may approximate Bonnor- 
Ebert structures which are however of variable mass and may or may not experience gravitational 
collapse, in qualitative agreement with the large observed frequency of cores with BE- like profiles. In 
our collapsing simulations, a time 0.5 Myr typically elapses between the formation of the core and 
the time at which it becomes gravitationally unstable, and another ~ 0.5 Myr are necessary for it to 
complete the collapse. 

Subject headings: ISM: clouds — ISM: evolution — ISM: structure — stars: formation — turbulence 



1. INTRODUCCION 

The process by which a gas parcel ("core") within 
a molecular cloud (MC) initiates a collapse leading 
to the formation of a star or group of stars remains 
loosely understood, in particular the details of its dy- 
namical evolution. Observations indicate that "prestel- 
lar" molecular cloud cores (i.e., those that do not 
yet contain a protostellar object, but that appear to 
be on route to forming it) have a dens ity structure 
that resembles Bonnor-Ebert (BE) profiles (|Ebert II1955I : 
[Bonnor 1956), being nearly flat in their central re- 
gions, while approaching the singular isothermal sphere 
(SIS) profile n(r) oc at large radii. "Stellar" cores 
(those already containing a Class or Class I proto- 
stellar object), on the other hand, appear to have den- 
sity profiles c loser to that of t he SIS througho ut their 
volu me (e.g.. lAlves. Lada fc Lada l |200lt ICaselli et al. I 
120021 : iKirk et al. 1 120051 : iLee et al. I boOTl ) (see also the 



reviews bv iLada et al. 



I Ward-Thompson et al. 



120071 : Idi Francesco et all 120071 : 
12007 1. and references therein). 



The line profiles and spatial distribution of molecular- 
line observations provide further clues to the dynamics. 



For example, based on obser vations of C S (3-2) CS(2-1), 
DC0+(2-l) and N2H+(l-0), EeeetaD (l200l found a 
moderate fraction of prestellar cores (18 out of 70 in their 
Table 2) showing clear evidence of subsonic inward radial 
motions, at velocities v < 0.07km s~^. Moreover, studies 
of individual starless cores have suggested that the radial 
veloc ity does n ot increase ap preciably towards the center 
(.Tafalla et aTlfipoa: .Williams et al. |[i99a: Tafalla et al. I 
120041: ILee etH. Il2007l:"1s'chnee et al.ll2007F . Those inward 
motions frequently extend to long enough distances from 
the cores' centers (a few tenths of a parsec) that they 
se em inconsis tent with the "inside-out" collapse model 
of lSbll (fl977l ). since a central protostar should have had 
time to form b y the time the rarefactio n wave reaches 
those distances (|di Francesco et aO2007l and references 
therein). 

A large number of theoretical studies have inves- 
tigated the collapse process starting from a variety 
of initial and boundary conditions, both analytically, 
through si milarity solutions, and numeric ally (e.g. , 
Larson 1196 9: Penston 1969 allbl:[Shuiri 977: Hi mter Ill977t 



Foster fc Chevalier ill993t iHennebell e et al. 1120031 ). All 
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of these studies have considered the collapse of a fixed 
m^ass of gas, either through the usage of a hot, ten- 
uous confining medium that pressure-confines the core 
while adding no weight to it, or through fixed bound- 
aries. Moreover, most of these studies used static initial 
configurations, either with uniform density or with BE 
hydrostatic equilibrium profiles. 

On the other hand, MCs are thought to be su- 
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personicall y turbulent, since they ex hibit supersonic 
linewidths (jZuckerman k. Palmer iflQT^ ). and MC cores, 
as well as their parent MCs themselves, have been 
suggested to be turb ulent density f luctuations within 
their envir onmen ts (|Sasao ' 1973; lElme green 
iBallesteros-Paredes et al. I [1999) . bei ng produced by 
effecti vely supersonic compressions. iHunter fc Fleck I 
(|1982| ) showed that the effective Jeans mass of a fluid 
parcel subject to an external compressive velocity field is 
significantly decreased with respect to its norma l static 
value. Furthermore, IVazquez-Semadeni et al. I (|2005D 
have recently pointed out that, if MCs are isothermal 
throughout,'* then the hot, tenuous medium necessary to 
confine and stabilize a hydrostatic equilibrium configu- 
ration is not available, and the equilibrium state is then 
expected to be unstable in general. 

This leads naturally to the question of whether hydro- 
static equilibrium configurations can be produced in such 
turbulent conditions, and if so, how do they arrive at that 
state. Otherwise, if the entire process is dynamic, one 
can ask what is the density and velocity structure of the 
cores at the time they engage into collapse and, if differ- 
ent from the initial conditions normally assumed, what 
effects does it have on the evolution. Moreover, if the 
core is formed and induced to collapse by a compressive 
wave, then in general there is an inflow that builds up the 
core dynamically, and the mass that ends up collapsing is 
not previously determined by the initial conditions, but 
rather is determined "on the spot" depending on the local 
instantaneous conditions. The studies of collapse men- 
tioned above cannot answer these questions, since they 
already assume gravitationally unstable structures, and 
initial hydrostatic equilibrium configurations, so that all 
of the mass is involved in the collapse. 

A study that comes close to these goals is that by 
iHennebelle et al."1 ()2003l ). who investigated the effect of 
increasing in the pressure external Poxt to an initially 
stable BE sphere. They noted that the resulting config- 
urations are a good match to the observations because 
the density profile is flat at the center, and the prestellar 
phase is characterized by subsonic inwards velocities at 
the outskirts, and by nearly zero velocity at the inner 
parts. However, having a hot conflning medium outside 
and an initial hydrostatic proflle, this study still could 
not capture the core formation part of the evolution, 
and predetermined the mass that collapses from the ini- 
tial conditions. Also, it did not consider the possibility of 
a transient compression and thus of a failure to collapse. 

A brief discussion of t he dynamic scenario of co re for- 
mation has been given bv lWhitworth et al. I (|2007f ) in the 
context of the formation of cores that give birth to brown 
dwarfs. These authors have suggested that the dynamic 
formation of cores should involve a mass growth period 
and the conflnement of BE-like structures by ram pres- 
sure of external inf ailing material. 

In view of the above, in this paper we then present nu- 
merical hydrodynamical simulations in spherical coordi- 
nates of transient compressions in homogeneous, initially 
gravitationally stable regions, with the purpose of inves- 

* The possibility that MCs may contain warmer atomic gas 
interspersed with the colder molecular gas has been recently raised 
by Hcnnebelle & Inutsuka (2006). If this turns out to be the case, 
then the argument against the feasibility of stable hydrostatic BE- 
like MC cores weakens significantly. 



tigating the formation of cores embedded in turbulent 
molecular clouds. In particular, we focus on the evolu- 
tion of its density and velocity proflles, the timescales 
required for a core to be assembled and then collapse or 
redisperse, and the mechanism by which a certain frac- 
tion of the mass is gravitationally "captured" to then 
proceed to collapse. 

In particular, the timescale issue is highly relevant be- 
cause it is often thought that the prestellar lifetimes of 
the cores in the turbulent scenario of star formation are 
of the order of one core's free- fall time ts- H o wever , 
it h as been shown by [V azquez- Semadeni et"an (j2005l ) 
andl Galvan-Madrid et al. ■ (,2007|) that even in highly dy- 
namical, driven-turbulence simulations, the lifetimes are 
a few to several times Tg. It is important then to in- 
vestigate the detailed evolution of cores formed by tur- 
bulent compressions, to understand the reason for those 
observed timescales. 

The plan of the paper is as follows. In ij2] we first 
discuss the motivation and applicability of the spherical 
symmetry used in this paper, and then we describe the 
numerical setup of the problem. In fj3] we present the 
results of two fiducial cases of core evolution, one col- 
lapsing and one rebounding, and in SJJ] we then disccus 
the implications of our results, and compare with exist- 
ing observational and theoretical work. Finally, in ^ we 
present a summary and some concluding remarks. 

2. THE MODEL 

2.1. The need for focused compressions 

The formation, and subsequent induction to grav- 
itational collapse, of clouds and clumps by com- 
pressive velocity fields (as opposed, in particular, 
to collisions of pre-existing clouds) has been stud- 
ied intensively by numerous worke rs for more tha n 
three decades ( e.g., S eisao J973; Hunter & Fleck 1119.8 
Tohli ne et all 119871: TMcKee et al. I 1993; Elmcgrcen. 
19931: IVazauez-Semadeni et al. Ill996l ,2003. .2005j . .200 
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Klesscn et al. ll20'o a 'Heitsch et al. ll2001l:lLi et al. I 
Tilley & PudntzJ ,2004, 2005). In such scenario, 
the formation of cores and stars is in agreement with 
observational studies that suggest that s tar formation 
is a rapid and dynamic process ( e .g., iLee &: Mverd 
'1999- 'Ballesteros-Paredcs et al.' '1999'; 'Elmcgrcen' '20001 
Prin g lc et al. 2001; Briccno et al. 2001; Hartmann et al] 
l200lTiBallesteros-Paredes fc Hartmannll2007D . The abil- 
ity of a compression to induce collapse is directly related 
to the stability of self-gravitating equilibrium structures, 
which in turn depends critically on the geometry of the 
configurations and on the effective polytropic exponent 
(7off) of the medium. This exponent describes the re- 
sponse to compressions of a medium subject to heating 
and coohng processes (iTbhline et aTl 119871 : lElmegreenl 
Il991t IVazquez-Semadeni et al. 1119961 ). so that the flow 
exhibits an effective polytropic equation of state of the 
form P oc p'^"" . 

The stability of self-gravitating structures depends 
both on the geometry and on 7eff because both of them 
influence the variation of the ratio = \Eg\/Eth upon 
compressions, where \Eg\ is the absolute value of the 
gravitational energy and Eth is the (supporting) internal 
energy. For example, it is well known that the existence 
of stable spherical configurations without any external 
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confi ning agent requires 7eff > 4/3 ([Chandrasekhar I 
fl98ll §117). In this case, decreases upon a com- 
pression, and increases upon an expansion, rendering 
the equiUbrium stable. This behavior is reversed for 
7off < 4/3, so that equihbrium configurations are un- 
stable in this case. We refer to the value of -fcs at which 
the reversal occurs as the critical value, 7cff,c- 

Now, if the compression occurs along ly directions, so 
that the density increases as L", where L is the length 
scale along the direction (s) of compression, then the 
rate of variation of |i?g| and of with the compres- 
sion depends on v and, as a consequence, the critical 
value of 7eff also depends on v. Specifically, on e obtains 
(jMcKee et al. II1993I: FVazguez-Semadeni et al. II1996D 

7cff.c = 2(l-i.-i). 

This recovers the value 7cff,c — 4/3 for three-dimensional 
(e.g., spherically symmetric) compressions, as well as the 
well known fact that planar compressions {v = 1) can- 
not induce collapse in isothermal flows (7off = 1), since 
7eff,c = in this case. Equation ([1]) also shows that the 
induction of collapse in isothermal media requires com- 
pressions in more than two dimensions {v > 2). 

These considerations show that useful insight can be 
gained from the analysis of spherically symmetric com- 
pressions, like those assumed in this paper, since the 
compressions that induce the collapse of selected sub- 
regions of turbulent, isothermal molecular clouds need 
to be of dimensionality higher than 2. Such compres- 
sions are expected to be rare but still existing in general 
supersonic turbulent regimes, involving a certain degree 
of focusing (or converge nce) of the flow. For example, 
IWhitworth et al. I (|2007D appeal to the low but flnite 
probability of such focused compressions to explain the 
scarcity of ^ O.OIM0 brown dwarfs. 

Finally, note that these convergent flows are in princi- 
ple different from a simple passing shock (which is essen- 
tially a planar compression), although the latter can also 
induce multidimensional compression when the shock has 
a finite transverse extent, producing a flattened structure 
that then can contract gravitationally in the tranverse 
direction. 

2.2. The numerical setup 

In view of the above considerations, in this work we 
consider an idealized spherical cloud subject to an ex- 
ternal compression wave. The simple spherical geometry 
allows us to focus on the basic phenomena related to the 
effects of the compression. Moreover, by considering a 
uniform density distribution and allowing the system to 
dynamically choose the amount of mass involved in the 
collapse, we avoid some of the restrictions that previous 
work has imposed on the evolution. 

The hydro dynamic evolution of th is setup was solved 
using ZEUS (jStone fc Normanl[l99l l. a finite difference, 
time explicit, operator split, hydrodynamic code. The 
calculations were performed on a ID spherical grid, with 
the domain spanning the range < r < 1 pc with 1000 
grid points spaced such that 5ri^i/5ri — 1.005. This 
yields a spatial resolution of « 3 x 10~^ pc at the in- 
ner boundary, and « 5 x 10^'^ pc at the outer bound- 
ary. (Selected simulations were performed with a much 
higher resolution of 4000 grid points and no significant 
differences were observed.) The boundary conditions are 



"reflecting" at r = and "outflow" at r = 1 pc. No 
confining agents are used whatsoever (closed boundary 
nor hot tenuous medium) , implying that mass can freely 
leave the system, although it cannot enter. The absence 
of a confining agent attempts to emulate the situation 
of a density enhancement immersed in a much more ex- 
tended medium at the same temperature. 

All simulations started with a constant density distri- 
bution, an isothermal equation of state, and a temper- 
ature T = 11.4 °K which, with a mean particle mass 
/i = 2.36 mu, yields an isothermal sound speed Cs — 
0.2 km s~\ This setup was perturbed by a compressive 
velocity pulse given by the relation 

{0, r<ro-dro; 
-«osin(fi^), ro-rfro<r<ro + dro; 
-vo sin (1^^) , ro + rfj'o < r\ 

(2) 

where vq and rg are parameters of the simulation, dr^ — 
O.lpc, ri = [rinax + (ro + d?-o)]/2, and dri = [r^^^-{ro + 
c?ro)]/2, with r-max — Ipc A simple self-gravity module 
was also added to the code. 

Our approach continues along the lines of simple, basic 
models that ha ve explored the gra vitational c ollapse of 
MC cores, since La rson I (1 19691) an d iPenston I ()1969 a"b^ , 
through [Hennebelle et al. I (|2003D . which we have ex- 
tended to include an initial velocity impulse, intended 
to mimic the random compressive motions expected in 
a turbulent medium. Nevertheless, the one-dimensional 
nature of the model, together with the adopted spher- 
ical geometry, makes this setup somewhat unphysical 
as it restricts the nature of the compressible wave to 
spherical shells, and "turbulent" support in this model 
is present only as purely divergent motion, with no ro- 
tational component. A more realistic way of modeling 
the core formation process, albeit perhaps less amenable 
to detailed analysis, would be to perform full 3D nu- 
merical simulations, via random compressions of finite 
cross-section generated by bulk motions of the gas, sim- 
ilarly to what has been d o ne for the diffuse medium by 
I Vazquez- Semadeni et al. I (|2007t ). We intend to pursue 
this in the near future, over the theoretical foundation 
laid out by the simple present study. 

Another limitation introduced by the adopted geome- 
try is the large mass of the collapsed core resulting from 
our simulations (cf. ij3.2p . In a more realistic simula- 
tion, without the geometrical and symmetry restrictions, 
the collapsing system would probably undergo fragmen- 
tation. Therefore, we see the collapsed objects generated 
in these simulations not as a single star, but as the pre- 
cursors of small clusters. 

3. THE SIMULATIONS 

3.1. Spontaneous collapse 

In order to study the effect of velocity fields in induc- 
ing the collapse of molecular cloud cores, we first need 
to determine when they can collapse under the influence 
of their self-gravity alone. Because of the adopted spher- 
ical geometry (the usual Jeans analysis is applicable to 
sinusoidal perturbations in plane parallel geometry) , the 
critical density pc and mass (which we refer to as the ef- 
fective Jeans mass) at which the core collapses may differ 
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Fig. 1. — Spontaneous collapse of the cloud without velocity impulse. Left: when started with a constant density n = 175.70 cm~'^, 
the central region of the cloud undergoes gravitational collapse after a small bounce off the center. Right: when started with a constant 
density n = 112.72 cm""^, the cloud bounces off its center and expands until it disperses. 



slightly from the standard Jeans values, and so we deter- 
mine them here numerically. We set vq = and let the 
simulation run for lOMyr with a series of different initial 
densities. 

As the simulations are started, self-gravity causes the 
cloud to begin contracting, increasing its mean density 
(see fig. [1]). At some point, the pressure gradient in the 
inner parts stops this process and the contraction is re- 
versed (the cloud "bounces" momentarily). If the cloud's 
mass is large enough, self-gravity takes over again, the 
expansion is also reversed and the cloud collapses; other- 
wise, the expansion continues until the simulation ends. 
It is found that an initial density value of 160 cm~'^ yields 
a collapsing core, while a 2% lower density does not; 
therefore, we take the critical density as pc — 160 cm~^. 
At this density, the mass in our numerical box (of ra- 
dius R — 1 pc) is Mbox(Pc) = 39.1M0. For comparison, 
the mean density for which the standard Jeans length 
equals the diameter of the numerical domain (2 pc) is 
p = pj= TTc^/GLj = 125cm-3. 

In the light of this result, we define the effective Jeans 
mass as the spherical Jeans mass (i.e., a sphere with di- 
ameter equal to the Jeans length at mean density p) times 
a fudge factor A so that the product equals the box's 
mass at the empirical critical density: 

where is the Boltzmann's constant, T is the temper- 
ature, fj, — 2.36 mu is mean particle mass, and G is the 
gravitational constant. By setting Afj.cff = Mbox, we ob- 
tain A — 1.45. For comparison, the standar d Jeans mass 
at Pc is Mj = 27.0Mq, and the BE mass (|Ebert I HQS! 

IBonnor ifTOSll is Mbe = l-lSc^/ {G^Pc)^^^ = IO.OMq. 
3.2. Cores formed by ram-pressure 



Although a large number of simulations were per- 
formed, our discusion will focus on two of them, re- 
spectively representative of non-collapsing and collapsing 
cases. We shall call SI the simulation with the initial ve- 
locity impulse at tq = 0.33 pc, while simulation S2 places 
the impulse at ro = 0.67 pc. Both simulations have the 
same velocity amplitude {vq — 0.4 km s^^ = 2cs) and 
sub-critical initial density (112. Tcm""^ w O.Tpc), meaning 
that in the absence of compressive motions, both simu- 
lations would simply expand away. 

The evolution of simulation SI is shown in Figure [H 
This figure respectively shows, as a function of radius, 
the density, the logarithmic slope of the density radial 
profile, the velocity, and the core's mass (solid line) and 
Jeans mass (dashed line) inside the radius, in the four 
rows from top to bottom at selected times (left to right 
columns). Shortly after the starting time (t w 0.28 Myr), 
a shell bounded by two shocks appears on the inner side 
of the initial velocity pulse, at logr « —0.7 (r w 0.2 
pc) . The formation of these two shocks and the shell be- 
tween them is due to the middle parts of the compressive 
wave, which have the highest velocities, catching up with 
the frontal parts of the wave, causing the formation of 
a shock, which splits into two shocks receding from each 
other and leaving the shock-bounded layer in between. 

The innermost one of these shocks propagates inward, 
leaving a large amount of inflowing mass behind it (fig. 
[2^)). As it travels toward the center of the cloud (t « 0.77 
Myr), geometrical focusing dramatically increases the in- 
ternal density, lowering the effective Jeans mass of the 
inner parts of the cloud (fig. [2]c) . As the shock bounces 
off the center and expands outwards, the shocked gas 
behind it is left at uniform density and at essentially 
zero velocity; that is, a quiescent core is formed, with 
the shock-bounded shell dubbing as the core's envelope 
(d). The gas from the shock-bounded shell continues 
to fall in, being incorporated into the quiescent core as 
it passes through the inner shock. Although the mass 
of the quiescent core increases, its density is somewhat 
lowered because of a mild expasion of the compressed re- 



5 




-4 -3 -2 -1 0-4 -3 -2 -1 0-4 -3 -2 -1 0-4 -3 -2 -1 0-4 -3 -2 -1 

log r [pc] 



Fig. 2. — Evolution of simulation SI. Each row shows the evolution (down from the top) of density (n), logarithmic den- 
sity slope (d log n/d log r), velocity (fr), mass internal to radius r {M{r), solid line) and effective Jeans mass (Afjcff, dashed line) 
at 0.000 (colum n g), 0.625 (b), 0.775 (c), 0.925 (d), and 1.500 Myr (e). Arrows show the direction of motion of the shocks. [See 
|http : //wHW. astrosmo ■unam.mx/'^g. gomez/publlca/f 2 .mpg for an mpeg animation of this figure.] 
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Evolution of simulation S2 at 0.000,1.125,1.525,2.000 and 2.625 Myr. 

|http : //www. astrosmo ■unam.mx/'^g. gomez/publica/f 3 .mpg for an mpeg animation of this figure.] 
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gion. As a result, the Jeans mass becomes larger in the 
innermost parts of the core, and decreases close to the 
shock-bounded layer and within it [compare panels (d) 
and (e) of fig. [2]. As the core acquires more mass, its 
density profile starts to deviate from being uniform, and 
to approach that of a truncated BE sphere. However, in 
this simulation the mass of the inner core never becomes 
equal to Mj^cff at any radius, and the uniform-density 
core begins to expand indefinitely, developing positive 
velocities at the outermost regions first (fig. 2e). 

In simulation S2, the early evolution is quite similar to 
that of SI. At i « 1.5 Myr, the inner shock bounding 
the layer bounces off the center of the core (fig. [Sjc) and 
begins traveling outwards. But then, some 0.5 Myr later 



{t w 2.0 Myr), the amount of mass in the core finally be- 
comes equal to Mj cff (rf) at r « 0.07 pc, and from that 
moment on, collapse ensues, culminating with the forma- 
tion of a singularity at i « 2.6 Myr (e). It is interesting 
that at i « 2.0 Myr, the mean density in the quiescent 
core is n « 2.85 X lO^cm"^, implying a free-fall time 
Tff « 0.2 Myr, less than half the time the actual collapse 
takes. There are several reasons possibly responsible for 
the discrepancy with the observed collapse time of 0.6 
Myr (from t = 2 Myr t o ^ = 2.6 Myr) . For example, 
it was already noted by iLarson I (|1969l . Appendix C) 
that the actual collapse time lasts nearly 1.5 times the 
free-fall time, because the pressure gradient is never neg- 
ligible. The remaining difference is probably due to the 
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imprecisions introduced by considering the mean density 
rather than the detailed radial distribution, and the fact 
that the core is increasing its mass, so that the instability 
sets in at an undetermined radius. 

A very interesting feature of both simulations is the 
fact that the density structure of the core+envelope sys- 
tem resembles a BE sphere during the period over which 
the shock front is traveling outwards (figs. [2]e and[3li), 
at which times the innermost parts of the core have a 
nearly constant density and the shock-bounded layer ap- 
proaches an r"^ density profile. After the formation of 
the singularity at the center, the density profile ex- 
tends throughout the core, similarly to an SIS profile (fig. 

Els). 

It is worth remarking that, even though the idealized 
geometry and initial conditions adopted in this paper 
should not have a strong impact on the general qual- 
itative behavior of these simulations, the above quoted 
sizes and time-scales are expected to depend on these de- 
tails. In fact, the outcome of the simulation (collapse or 
re-expansion) already depends sensitively on the param- 
eters of our idealized simulations. For example, the only 
difference between simulations SI and S2 is the initial 
position of the compressive pulse. Moreover, a simula- 
tion similar to SI but started with vo — 3cs and a lower 
initial density (88.8 cm~^) collapses 0.4 Myr after the ini- 
tial shock bounces off the center. That is, a lower initial 
mean density can be counterbalanced by a larger com- 
pressive Mach number, in the quest for inducing collapse. 
Nevertheless, this simulation still goes through the same 
qualitative evolution as simulation S2. For these reasons, 
the results presented in this section can be regarded as 
a qualitative description of the formation and evolution 
of molecular cloud cores, while more quantitative anal- 
ysis requires more realistic simulations involving multi- 
dimensional compressions, and it is left for future work. 

4. DISCUSSION 

4.1. Generality of shock-bounded self- gravitating 
structures 

The formation of growing shock-bounded structures is 
not exclusive to the spherical symmetry used in this pa- 
per. Planar compressions are generally known to produce 
shock-bounded layers in both isothermal and radiatively- 
cooling flows. The one-dimensional plane-parallel prob- 
lem is equivalent to that of a shock front hitting a wall 
and then reflecting off it, and thus by construction the 
gas between the wall and the shock is at rest, with the 
shock front receding from the wall at the post-shock ve- 
locity of the flow. The shocked layer increases its mass 
as gas from the incoming flow is incorporated into it af- 
ter crossing the bounding shock. The shock bounded 
layer is thus the planar equivalent of our spherical-shock- 
bounded quiescent core. The plane-parallel problem has 
been worked out analytically in ID and numerically in 
two or three dimensions by Folini & Walder (2006) for 
the i sothermal case and bv | Henne bellc & Perault I (jl999l) 
and IVazquez-Semadeni et al. I " 12006.) for the thermally 
bistable case. The main difference with the spherical case 
is that plane-parallel compressions in thermally bistable 
media (s uch as the warm HI gas) can induce gravitational 
collap se (|Hunter at al. 1 1 19861 IVazquez-Semadeni et al. I 
|2007() because in this case the flow behaves effectively 
as if having jcS < 0, while in an isothermal medium 



focused (i.e., multidimensional) compressions are nec- 
essary because in this case jcS — 1 (see §2.ip . How- 
ever, the planar and the spherically symmetric cases 
are qualitatively similar in that both involve the forma- 
tion of a shock-bounded structure that grows in mass 
by accretion through the shock until it becomes gravi- 
tationally unstable and begins contracting. In the pla- 
nar compression case, this process h as been modeled nu- 
merically by various workers ( e.g., iHunter at al. I fl986l : 
IVazquez-Semadeni et al. Il2007[ ). 

4.2. Implications 

The evolution of simulation S2 has a number of inter- 
esting important implications, which we now discuss. 

First, a compressive wave (or a negative-divergence ve- 
locity field) does not directly induce the collapse of an 
initially sub- Jeans core. The collapse happens only if 
at some point in the evolution the mass becomes larger 
than the Jeans mass. In all the collapsing simulations we 
have performed here, this occurs only after the resulting 
shock front has rebounded off the center, traveled out- 
wards, and incorporated a large enough amount of mass 
into the central core, so that a "traditional" Jeans crite- 
rion for collapse [M(r) > Mj^off] is satisfied there. Since 
the material behind the shock is left at zero velocity, no 
turbulent support is ever at play there. That is, the col- 
lapse does not occur because turbulence is dissipated in 
the core, as it is often believed, but rather because the 
growing core eventually reaches the effective Jeans mass. 
Moreover, as the shock continues to move outwards, the 
size of the region acquiring the effective Jeans mass in- 
creases, so that the determination of the mass that is 
subsequently incorporated into the collapse happens "on 
the spot" in a highly fortuitous manner. 

Second, a near r~^ density profile is approached at 
late times in the infalling envelope around the central 
core, both in collapsing and non-collapsing cases. The 
central core, in turn, evolves from a near-flat density 
profile to that of a truncated BE sphere as its mass in- 
creases. The central core and the envelope are separated 
by the outwards-traveling shock, which is, however, very 
mild, with a Mach number very close to unity. Thus, the 
core-|-envelope system may easily be taken for a single 
structure with a density profile resembling that of a BE 
sphere. 

Third, the velocity profile of the central core at all 
times after it has formed, has nearly zero velocity 
throughout. This provides a physical basis for the ex- 
istence of quiescent (subsonic non-thermal velocity dis- 
persion) and coherent (non-thermal velocity dispersion 
nearly constant through the core) cores, which in the tur- 
bulent model of core formation are the st agnation points 
of th e turbulent fiow in molecular clouds (jKlessen et al. I 
[20051) . 

Fourth, there is a time delay between the formation 
of the core and its gravitational collapse. The quiescent 
core grows from the center as the shock moves outwards 
incorporating mass into the central shocked region. This 
process cannot happen instantaneously, but rather re- 
quires a finite time until the core's mass equals the ef- 
fective Jeans mass. In our simulations, roughly 1 Myr 
spans from the moment of central core formation to the 
development of a singularity at the center, with roughly 
half of it being spent without any tendency to collapse. 
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This time delay naturally explains the high frequency of 
observed starless cores with BE-like profiles. 

Fifth, it is important to remark that even though 
our quiesecent cores are morphologically similar to BE 
spheres, they are dynamically very different: they are 
not confined by the thermal pressure of a hot, tenuous 
medium, but instead are confined by the ram-pressure of 
the inflowing gas from the envelope, and growing in size 
and mass accordingly, until they become dominated by 
gravity, at which point they engage into collapse. 

Sixth, and finally, it appears that the whole evolution 
is not very amenable to a similarity solution because: a) 
the initial velocity pulse is finite, so that the external flow 
is not rescalable. b) The inner shock bounding the shock- 
bounded shell hits the center and bounces back towards 
the exterior so the time of collision at the center breaks 
the self-similarity, c) The central core gradually increases 
its self-gravity and eventually may become gravitation- 
ally unstable, a process that continously transforms the 
core's density profile from uniform to being BE-like, first 
stable and then unstable. Similarity solutions may be 
most applicable after the format ion of the ce ntral singu- 
larity, as originally suggested bv lShu I (|1977l ). 

4.3. Comparison with previous work 

It is interesting to put the results of our numerical sim- 
ulations in context with those of previous studies. The 
main difference is that our simulations have investigated 
the formation of the cores in addition to their subsequent 
collapse, in order to study whether BE-like structures can 
be spontaneously produced out of supersonic turbulent 
compressions in isothermal molecular clouds. Thus, in 
particular, our study sheds light on the realizability of 
the initial conditions used by previous works. 

Our results suggest that in fully isothermal molec- 
ular clouds (i.e., without a warm, tenuous interclump 
medium that can stabilize a density enhancement), col- 
lapsing structures formed by random turbulent com- 
pressions in the medium morphologically resemble BE 
spheres through a large fraction of their evolution, be- 
cause they consist of a central core and an infalling enve- 
lope which, at late times after the formation of the cen- 
tral core, has a density profile with a slope close to —2. 
This extends previous results that plane-of-sky angular 
averaging and line-of-sight averaging cause the observed 
density profiles to be smooth er that the actual ones and 
easily confused with BE ones (|Ballesteros-Paredes et al71 
l2003t [Hartmannll2004 ) . Also at late times, near the onset 
of gravitational collapse of the core, the latter also devel- 
ops a density profile close to that of a BE sphere, which 
connects with that of the envelope. This means that, 
at the onset of gravitational collapse, our simulations fa- 
vor initial conditions for c ollapse such as those used by 
[Foster fc Chevalier I (|1993D . albeit with the added ingre- 
dient of a continuous accretion at the bounding shock. 

The establishment of a near r~^ density profile in the 
envelope at l ate times is i nteresting in the context of the 
discussion bv lShu I (|1977| ). He points out that the devel- 
opment of such a profile requires that the initial motions 
in the outer regions of collapsing cores be subsonic, so 
that all fiuid parcels are in acoustic contact with each 
other, and can therefore approach detailed mechanical 
balance. In our simulations, however, this need for an ini- 
tially subsonic condition appears to be in contradiction 



with the supersonic nature of the initial pulse. However, 
the acoustic contact is restored in the envelope because 
it consists of shocked gas that has been thermalized and 
thus initially subsonic. The fact that a near-r~^ profile 
develops in the envelope even in the non-collapsing sim- 
ulation can be understood because the compressive pulse 
effectively removes the support for the outer layers, anal- 
ogously to the effect of an inside-out collapse. 

Some authors have already studied spherically sym- 
metric flows with shocks in the context of protostar 
formation using similar i ty methods (e . g.. | Shen &: Lou I 
l20f)llLou fcGaxni2(M [Lou fc Wangll2006[ ). Similaritv 
studies are extremely useful in extracting the underlying 
asymptotic behavior of real flows. Therefore, it is impor- 
tant to compare our numerical solutions with existing 
similarity solutions of self-gravitatin g clouds in th e pres- 
ence of shocks, in particular those of lShen fc Lou I (|2004 
herefater SL04), whose study most resembles our numer- 
ical setup. These authors presented two possible classes 
of self-similar shocked flow in the context of the dynam- 
ical evolution of protostars, depending on the asymp- 
totic behavior of the solutions near the center of cloud. 
Their Class I solutions had negative (inflow) velocities 
(cx — r^/^), a density profile p oc r"'^/^, and finite mass as 
asymptotic limits at r ^ 0, while their Class II had pos- 
itive (expansion) velocity (oc r), constant finite density 
and vanishing mass (oc r^) as the asymptotic behavior in 
the same limit. In both classes, an outward-moving shock 
separates a collapsing (or expanding) inner part and an 
accreting outer part. None of these behaviors are seen in 
our simulations at any time. Their Class II is similar to 
our solutions during the core-growth stage, in that it has 
a uniform central density and an accreting outer part, 
which has a counterpart in the infalling shock-bounded 
layer in our models. However, in our system the central 
core is neither expanding nor contracting, but rather it 
is at rest. This difference is most likely a consequence 
of self-gravity being neglegible in our cores during the 
early stages of their evolution. That is, unlike the SL04 
solution, where self-gravity is important at all radii and 
at all times, in our simulations the relative importance of 
self-gravity increases secularly with time, going from be- 
ing zero at the time of core formation to being dominant 
at the time when gravitational instability sets in. 

Another rece nt study that is closely related to ours is 
that by Henncb elle et al.1 (j2003l ). who numerically inves- 
tigated the effect of increasing the pressure external Pcxt 
to an initially stable BE sphere. These authors found 
that slow rates of increase of Pext cause the sphere to 
approach instability quasi-statically, but higher rates of 
increase produced a compressive wave that triggers an 
outside-in collapse. It is noteworthy, however, that they 
do not report the bounce of the compressive wave from 
the center that we find. This is most probably because, 
in their case, the wave compresses a previously-existing 
core that is in a (fragile) stable hydrostatic equilibrium 
state, and so the role of the wave is to directly trigger 
the collapse. Instead, in our case, the compressive wave 
forms the core, and adds mass to it until it becomes 
gravitationally unstable and proceeds to collapse. More- 
over, in the case of Henncbelle ct al. (2003), the mass 
of the core was fixed, being bounded by a hot, tenu- 
ous medium, while in our case, the fraction of the mass 
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that is driven to collapse is determined "in real time" by 
the interplay between the accreting gas and the outgoing 
shock wave, and moreover the mass that becomes grav- 
itationally unstable increases with time, so the collapse 
proceeds "inside out" but over an intermediate range of 
radii. Thus, we see that the choice of equilibrium or 
out-of-equilibrium initial conditions and continuous or 
discontinuous boundary conditions leads to very differ- 
ent patterns of evolution. Which model applies best to 
actual turbulent molecular clouds probably depends on 
whether they consist of a single, nearly isothermal molec- 
ular phase (our model), or of a mixture of colder, denser 
molecular cloudlets im mersed in a more tenuous an d 
warmer atomic medium (jHennebelle &: Inutsuka II2OO6D . 
Extensive theoretical and observational work, focusing 
especially on the velocity structure of the core s, is needed 
to dec ide on this issue. The recent results of iLee et al. I 
(|2007f ). indicating the presence of a sharp infall veloc- 
ity increase at ^ 0.03 pc from the centers of the starless 
cores L694-2 and L1197, would seem to favor our dynam- 
ical scenario for the formation of the cores. 

Finally, our res ults are fully cons is tent with the sce- 
nario outhned bv IWhitworth et al. I (j2007l ). These au- 
thors have foreseen the formation of evolving BE spheres 
bounded and fed by the accretion of external infalling 
material, which can collapse if the core eventually reaches 
the BE mass. Although they restricted their discussion 
to the formation of brown dwarf-producing cores, the 
simulations described here are seen to be applicable to 
the formation of cores of arbitrary mass. It is indeed 
likely that, as the core to be formed is of smaller mass, the 
required compression and focusing need to be stronger, 
as the initial conditions will have sizes much smaller than 
the local Jeans length (cf. §3.2p . 

5. SUMMARY AND CONCLUSIONS 

In this paper we have performed a numerical study 
of the formation of dense cores by dynamical compres- 
sions in isothermal, non-magnetized media, using simple 
one-dimensional calculations in spherical geometry. Our 
results show that cores assembled by this process consist 
of a central, quiescent core with density 10^ cm~'^ that 
grows in mass and size as it accretes mass from a sur- 
rounding envelope. The quiescent core and the envelope 
are separated by a mild shock with Mach number just 
above unity, and the accretion from the envelope provides 
ram-pressure that confines the central quiescent, growing 
core. As the central core increases its mass, it passes first 
through a neglegible-gravity, uniform-density stage, and 
later, as self-gravity becomes important, it evolves into a 
"pseudo BE-sphere" stage. If at some point in the evo- 
lution, the mass of the core-envelope becomes lager than 
the Jeans mass, the core proceeds to collapse. Otherwise, 
it begins to reexpand. Even in collapsing cases, this pro- 



cess requires a relatively long time to complete, taking 
'^0.5 Myr from the first appearance of the central core 
to the time at which it becomes gravitationally unsta- 
ble, and another ^ 0.5 Myr for the collapse to produce 
a singularity at the center. 

At all times after the formation of the central core, 
the combined density structure of the core-|-envelope 
system resembles that of a BE sphere, since it is flat- 
tened at the center, and has an envelope that approaches 
an density profile at late times. Moreover, the 
central core is quiescent at all times, except for the 
very late stages of the non-collapsing case, in which the 
core begins to expand. Thu s, the high observed fre- 
quency of BE-like profiles (e.g. Idi Francesco et al"~ll2007t 
iLada et al. 1 12007[ ) and quiescen t /coherent veloci t y dis- 
persion struc ture (iM versI Il983t iGoodman et al.l Il998t 
Caselh et alTI [2 002: T afalla et all 120021 : iTafalla et al. I 



20041: TSchnee et al.i .20071 ) is naturallv accomodated in 



this scenario of dynamic assembly of MC cores, as sug- 
gested also by studies of dense cores in turbulent simula- 
tions of 3D, isothermal molecular clouds .Klessen et all 
((2()05). However, the structures are not classical BE 
spheres, because they are confined by ram-pressure, 
rather than by thermal pressure, and are consequently 
accreting mass and growing in mass, size, and self- 
gravitating energy, in a process qualitatively similar 
to that described for th e formation of giant MCs by 
IVazguez-Semadeni et al. I ()2007[ ). In both cases, there 
is a secular evolution, characterized by the mass increase 
of the cloud or core. 

The velocity structure of the cores formed in our simu- 
lations appears consistent with recent radiative transfer 
models fo r the structure of cores L694-2 and L1197 pre- 
sented bv lLee et al. I ()2007| ). which exhibit a nearly zero 
central velocity and a sharp rise at radii ^ 0.03 pc. We 
plan to carry out a radiative transfer study of the den- 
sity and velocity structures produced by our models in 
the near future, in order to perform detailed comparisons 
with ob servational studi es based on multi-tracer stud- 
ies (e.g. iLee et al~ll2004l ) as well as on line-profile map- 
ping of prestellar core s (e.g., Tafalla et al. 1998. .200l 
ILee et al. lfT99l [20071: iSchnee et al.ll2007D . 
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